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I. INTRODUCTION 

Thermodynamic potentials such as the free energy are usually calculated in the infinite size 
limit, and in this limit only the leading (macroscopic) contribution is taken, as in studying 
bulk or macroscopic properties this gives a sufficient accuracy. Recent developments both in 
field theory and solid state physics have posed series of problems, where surface or impurity 
contributions are important [1,2]. As this contributions are less than the macroscopic ones, 
the otherwise formidable problem of calculating the next to leading order terms is set in 
focus. 

For the one dimensional (ID) systems soluble by Bethe Ansatz (BA) the free energy is 
calculated following the method developed by C.N. Yang and CP. Yang [3] for the 5 Bose 
gas. The basic idea of this method is that through the density of the momenta (rapidities) 
an entropy can be defined enabling one to write up the free energy functional of the finite 
temperature system. Minimization of this functional with respect to the momentum density 
yields the most probable distribution of the momenta and the value of the (macroscopic) 
free energy. In terms of the partition function, the minimizing of the free energy functional 
corresponds to finding the state entering with highest weight into the partition function 
(saddle point). The contribution of the states around this one (saddle-point fluctuations) 
is an 0(1) factor to the partition function, i.e. an 0(1) correction in the free energy (if 
not zero). Our aim is to develop a convergent functional integral method to calculate this 
contribution, what is actually a finite size correction to the free energy. 

First, to avoid special difficulties due to the BA, we write up, evaluate and discuss the 
convergence problems associated to the functional integral for the case of the free Fermi gas. 
Using the repulsive 5 Bose gas - the simplest BA system - as an example, we generalize 
the procedure to the BA systems. We give also a general expression for the contribution 
of the saddle point fluctuations in more complicated systems, and discuss the case of XXZ 
Heisemberg chain in more detail. 

The quantity we calculate for the free Fermi and 5 Bose gas is the grand canonical parti- 
tion function and the thermodynamical potential defined by its logarithm (grand canonical 
potential). In a strict sense this latter is different from the free energy, nevertheless it is 
not unusual to name it so. Throughout the paper also we use this certainly less precise, but 
hopefully not confusing name. 

We consider systems at periodic boundary conditions (PBC). As in this case there is no 
boundary, the naive conception, which identifies the 0(1) corrections as boundary con- 
tributions [1] would lead to the conclusion, that for PBC there are no 0(1) corrections. 
Surprisingly we find this is true only for the case of free Fermi gas, but does not hold for the 
interacting systems. In the case of the repulsive 8 Bose gas we calculate the contribution 
of the saddle-point fluctuations also for open end systems. We find that also in this case 
the corrections are given by bulk properties (the energy of the particles as a function of the 
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rapidities and by the kernel, i.e. the derivative of the scattering phase shifts), and the nature 
of the boundary is reflected in the structure of the contributions. 

We have to note, that our method works in the thermodynamic (infinite size) limit only, 
and as a consequence the correction we find gives the so-called residual entropy in the 
zero temperature limit [1] . For finite systems the entropy S tends to the logarithm of the 
ground-state degeneracy \n.D as the temperature T tends to zero. The residual entropy we 
find behaves differently: its T — > limit is not connected to the ground-state degeneracy 
D rather it depends on other properties like particle density or magnetization. This is a 
manifestation of the fact, that the size — > oo and T — > limits do not commute. 

The usual structure of the Bethe Ansatz equations makes it possible to give a general 
expression for the contribution of the saddle point fluctuations in a large class of Bethe 
Ansatz systems. Using the XXZ chain as an example we discuss some points one has to 
pay attention to when applying our formula. This considerations show that in addition 
to the saddle point fluctuations other effects can give 0(1) corrections too, more over in 
the isotropic (SU(2) symmetric) case the degeneracy of the states belonging to the same 
spin-multiplet can give an even larger contribution. 

The paper is organized as follows. In Sec. II we develop the method and discuss the con- 
vergence problems. We use the free Fermi gas to have a possibility to check the result. The 
method is generalized for the 8 Bose gas in Sec. III. Here we treat both the PBC case and 
the case of open ends with surface potentials. Sec. IV is devoted to the generalization of 
the method for more complicated systems. Our results are summarized in Sec. V. In the 
bulk of the paper we concentrate on the main line of our method, and the technicalities are 
collected in appendices. 



where L is the size and at PBC the J quantum numbers are integers. The energy of a state 
is 



II. THE FREE FERMI GAS 



The state of a free Fermion system is characterized by a set of wavenumbers 




(2.1) 



E = ^2 e(fci) where 




(2.2) 



with jj, being the chemical potential, by which the required particle number N is fixed [4], 
and the grand canonical partition function is 
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Z=]Texp }-pJ2e(h) . (2.3) 

{hi} I i J 

Here the summation over {ki} means summation over all possible ki sets, and (5 is the inverse 
temperature 1/T. 

Now we follow the method developed by C.N. Yang and CP. Yang [3] for the BA systems. 
We split up the k axis into intervals Ak, and introduce the density of ks (p(k)) and holes 
(ph(k)) so that the number of the wavenumbers and holes in the (k, k + Ak) interval is given 
by LAkp(k) and LAkp h (k) respectively. Obviously 

p(k) + p h (k) = ^- . (2.4) 
The number of states characterized by the same p{k) function is 

n[p(k)] = Uu( P (k)), with u{p(k)) = [ LM{ ^^ {k)) ), ( 2 - 5 ) 

where in Y\ k k labels the Ak intervals. In the following \nu(p(k)) is calculated by Stirling's 
formula 

= LA M PW) + + O ( T ^ m , TK L^ (2.6) 

where 

s(p(k)) = {(p(k) + p h {k)) \n(p(k) + p h (k)) - p(k) \np(k) - p h (k) \np h (k)} , (2.7) 



and 



<m) = -\ H ^ k) + \^E^±^. ( , 8) 



By means of uj(p{k)) a free energy functional is defined 



by which 



F[p(k)]=J2(Le(k)p(k)Ak-T\nuj(p(k))) , (2.9) 

k 



Z= ]T e- pF W<)\. (2.10) 

{P(k)} 
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(Here the summation extends over all possible p(k) distributions.) In the usual procedure 
\nu(p(k)) is calculated up to leading order in L only, i.e. 

\nu(p(k)) LAks(p(k)) (2.11) 

is taken, and the free energy functional 

F L [p(k)\ = L E (e(k)p(k) - Ts(p(k))) Ak (2.12) 
k 



is minimized with respect to p{k) leading to 

e- /Je[k) (2.13) 



Po(k) „-/3e(fe) 



Ph,o( k ) 
and 

F min = ~ E ^ (l + e-MQ) Ak = /~ In (l + e -*<*>) dk . (2.14) 



Now we can define the functional integral for Z. We evaluate the summation over the 
possible p{k) distributions through integrals over the p(k)s: as LAkp(k) (being the number 
of ks in the interval Ak) is integer, the correct approximation for the summs is 

E =► / ■■■ [U(LAkdp(k)) . (2.15) 

{ P (k)} k 

Introducing 

r(k) = p(k) - p (k) (2.16) 

we can expand F [p{k)] around po{k). In order to have the functional integral convergent, in 
u(p(k)) of (2.6) we have to take into account the terms next to the leading ones too (<^(p)). 
This gives 

Z = e -/3F min>< ( 217 ) 

Here in the exponent the first term comes from the expansion of Fl [po + r] , and the second 
term is <;(po)- Introducing new variables 



^ 2 po{k)p hfi {k) 
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one arrives at 



Ze f3F min = l[Q-J e -«W) 2 ^(A;)) = 1 , i.e. Z = . (2.19) 



We can build up confidence in this calculation, as it reproduces the exact result. As, however, 
we want to generalize the method to cases where no exact results exist, it is worth to examine 
it in detail: it involves several approximations and also an implicit cut-off procedure, and 
we should see the conditions under which these are valid. Actually we have to check two 
types of approximations: 

1. The application of the free energy functional. The question to be answered is whether 
the definition of the free energy functional is accurate enough. The partition functions 
(2.3) and (2.10) are certainly equivalent in leading order, but are not exactly equal, as 
(2.10) is obtained from (2.3) by a kind of averaging: in each Ak interval we take 

5>xpj -Pj2 e ( k j) | — > ex P {-Pe(k)Lp(k)Ak + \nu(p(k))} , (2.20) 

with J2{kj} meaning summation over all possible choices of Lp{k)Ak kjS out of the 
L(p(k) + p h (k))Ak possibilities, and k being a k value within the interval Ak. As 
we are interested in calculating corrections to the macroscopic free energy we have to 
define this averaging more precisely. In Appendix A we show, that taking for k the 
mean value of the ks in Ak we introduce an O ((Ak) 2 ) error to the free energy density, 
that disappears when the J2(- ■ -)Ak — > f(. . .)dk limit is taken. In a strict sense the 
mean value of the ks in Ak still depends on p(k) and pt(k), but in the practice one 
may take the middle of AA; as k. Although this introduces an 0(1/ L) uncertainty in 
the mean value of energy/particle in each AA; interval, as however, these are random, 
do not sum up to an 0(1) contribution. 

2. The approximations applied while evaluating (2.10): 

• Conditions for the application of Stirling 's formula. Each £ integral collects the 
main part of the contribution from a region —X < £ < X where X is a number 
of the order of 5-6 ( [5]). We may use Stirling's formula in the above form, if 
in the corresponding p region both LAkp(k), LAkp h (k) 3> 1. These lead to the 
requirements 



LAkp (k) \l-X 



2p h ,o(k) 1 
po(k) + ph,o(k) LAkpo(k) 



> 1 



(2.21) 



and 
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LAkp h , Q (k) \1-X 

i.e. both condition are satisfied if 

LAkpo(k) > 1 and LAkp hfi (k) > 1 . (2.23) 

As for large enough k po(k) — ► 0, one has to introduce a cutoff (say A) in the k 
space so that (2.23) are met for all \k\ < A, i.e. 

LA&po(A)>l. (2.24) 

This imposes a relation on the L and A. Similarly, if the temperature T is small 
enough Ph,o(0) — ► thus the method is applicable at large enough size but not 
at exactly zero temperature: the larger the size of the system is, the nearer the 
T = can be approached. 

We have to note, however, that since the entropy part of the macroscopic free 
energy functional is already obtained through Stirling's formula, the problem 
connected to the application of this approximation is in principle present in the 
calculation of the macroscopic part of the free energy too (Appendix D). 

• Limits of the £ integrals. The estimations (2.21-2.22) show also, that in case of 
(2.23) the limits of the £ integrals can be taken to ±oo. 

• The ?(po)- We have taken q(p) at po and we neglected ?'(po) r + s"(po)r 2 /2- The 
term linear in r shifts the center of the £ integral, while the quadratic term 
modifies the coefficient of £ 2 . It is not hard to see that these corrections in the 
exponent are O (l/LAkp (k), l/LAkp h:0 (k)) ones, what we may neglect if (2.23) 
holds. 

Now we may conclude, that our functional integral method is established in a strict math- 
ematical sense in the L — > oo limit only, and involves a cutoff procedure in the momentum 
space. We discuss the relation of the cutoff to the size and some details of the cutoff proce- 
dure (paying special attention to its connection to the macroscopic free energy) in Appendix 
D. 

This method is not applicable for T = 0, but the T — > limit is meaningful. (In the present 
case the T — > limit reproduces the T = result, but this is not so in the case of interacting 
systems, as we shall see later.) 

III. THE REPULSIVE S BOSE GAS 



2pp(k) 1 

Po{k) + ph,o(k) LAkp hfi (k) 



»1, 



^2.22) 



The case of periodic boundary conditions 
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The Hamiltonian defining the ID repulsive 5 Bose gas is 

dx 



^=-E7^ + 2c E<K^-^)> c>0. (3.1) 

i<j 



Using the Bethe Ansatz Lieb and Liniger has shown, that the diagonalization of this Hamil- 
tonian in a box with a length of period L can be reduced to the solution of the system of 
algebraic equations 

N fh. _ h \ N - 1 

Lki = 2nJ t - J2 2 tan- 1 ^ J , J, = — ^— (mod 1) , (3.2) 

with N being the number of particles in the system [6]. The energy of the system is again 
given by (2.2), but now the wavenumbers are determined by the Eqs.(3.2) instead of Eq.(2.1). 
The finite temperature description [3] goes like we treated the free electron gas with the 
difference, that due to Eq.(3.2) Eq.(2.4) is replaced by 

p(k) + p h (k) = a + J2 K{k, k')p(k')Ak' 

k> 

/oo 
K(k,k')p(k')dk' (3.3) 
-oo 

with 

a = — and K(k,k') = - — ? . (3.4) 

2tt 1 ; 2tt c 2 + (k - k'f 1 ' 

One should note, that only those p(k) distributions are meaningful (physical), for which this 
equation yields Ph(k) > for all k. The free energy is again given by (2.12), but now it is 
constrained by (3.3). Its minimization leads to 

Po(k) _ o - /3e ( fc ) 



Ph,o(k) 



e~^ k > , (3.5) 



with the energy e(k) determined by the equation 

e(k) = e(k) K(k, k 1 ) In (l + e~^) Ak' 



/oo . 
K(k, k!) In (l + e 
-oo ^ 



(3.6) 



dk' 



Once e(k) is found, po(k) is given by the equation 

1 r°° 1 

= TT^ a + loo TT^ K{k > k >^ dk ' ' ( 3 - 7 ) 
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and the minimal free energy is 



-Vln(l + e-^«)AA; = --^- /°° In (l + e^W) dA; . (3.8) 



Fmm (32n^'~y-~ I"" (32 



To continue we expand the (2.12) around po(k) and ph,o(k): 

FlW4)] ,^_ TE£A 4(M|^_^|_^) . (3 . 9) 

T 2 \Po(k) +Pm( /c ) Po(k)/ 

Here the quantities 

r(k) = p(k) - p (k) and r h (k) = p h (k) - p hfi (k) (3.10) 

are constrained due to (3.3) by 

r(Jfe) + r h (ib) = Y,K(h k') r{k')Ak' . (3.11) 

k' 

This way the functional integral for the partition function reads 

Z = e~ pF ^ J ■■■ Jl[ (LAkdr(k)) x (3.12) 



{-?H 



1 ' ' ' ' " l ,Po(fc) + p M (*) PhM p (k)J 2 p oPhfi LAk2n) f ' 



what is to be evaluated under the constrain (3.11). After eliminating r h {k) we arrive at 

y LAk l_ ( (r(k)+r h (k)f _ rl(kl _ r^)\ = 
k 2\p (k)+p h , (k) p hfi {k) p {k)J 

-I £ LAM*) W ^£^ m ^ k")r{k")Ak" 



(3.13) 



+ ]T LAkr(k) K ^[K (k')Ak' 

k,k> Phfi(k) 



■W L Akr\k) P -^±^. 
2^ y ' p (k)p h , (k) 



Now, just as in the case of free Fermi gas we introduce according to Eq.(2.18). This 
leads to 
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E f (*0 ( S «,k - M k , >k ) (<? fc , )fc „ - M k ,, k „) ffl) , 



(3.14) 



where £(k) is that of (2.18), 5 k;k i is the Kronecker symbol, and 

1 



po{k)p hfi (k)Ak 
Ph,o( k )\ po(k) + p h ,o(k) 



K(k,k') > 



Po(k>) Phfi (k')Ak> 
\ Po(k') + p h , (k') 



(3.15) 



Finally, changing the integration variable to £ in (3.12) we have 



Z = / • • • / II f — ) exp I - E W (5 k , k , - M Kk ,) (8 k „, k , - M fe „, fe ,) Z(k?) 

J J k \ 71 ) \ k,k',k" 



(3.16) 



e -/3F min (det [ Sk>k , - Kw])- 1 



with 



Kk,k> = \Jphfi{k) ■ M k:k , ■ 1 



p (k)Ak 



K(k, k'), 



y/phflik') \ Po(k) + p h ,o( k ) ' \ Po(k') + p h ,o(k') 



Po(k')Ak> 



(3-17) 



The determinant in the above formula evaluated using the identities 

(det [I - K])- 1 = exp {-Trln(I - K)} = exp |E -TrK7 

1 n n 



(3.18) 



and taking the J2k^k — > / dk limit yields 

1 — e -f3F min +AS 



(3.19) 



where 



as = y -ic n 



(3.20) 



with 



/oo roo 
■■■ dh---dk n x 
-oo J —oo 

Po(ki) 



Po(h) +phfl{ki) 



K(h,k 2 ) 



po(k 



Po(k r . 



(3.21) 



po(k 2 ) + ph,o(k 2 ) po(k n ) + ph,o(k n ) 



K(k n ,k 1 ) 



/oo /-oo j I 

-oo ' ' ' L dkl '-- dK TT^ K(ku fc2) IT^ 



1 



1 _|_ e l3e(k 2 ) I _|_ g/3e(fc„) 



i^(/c n , fci) 
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(We have to note here, that, although our notation suggests so, for T 7^ AS is not a 
correction to the entropy, nevertheless we use this notation, as formally it comes from the 
density of states.) 

Now, similarly to the case of the free Fermi gas we should make a kind of "validity test". In 
this, in addition to the questions discussed there (application of the free energy functional, 
and approximations in the evaluation of the partition function) one hast to see also, that 
the J2 n * s convergent. 

1. The free energy functional. While in the case of the free Fermi gas the application 
of the free energy functional introduced to the free energy density an O ((Ak) 2 ) error 
only, in the present case we are faced to an apparently more serious problem, which 
originates from the fact, that the Jj quantum numbers in (3.2) are either integers or 
half-integers depending on the parity of the particle number. Changing all Jj quantum 
numbers from integer to half-odd-integers or vice versa shifts all hi by 7r/L leading to 
a shift also in the free energy. For general p(k)s this shift can be of 0(1) suggesting 
that the free energy for this system as a function of the p{k) is defined with an 0(1) 
accuracy only. Fortunately not this is the case. In Appendix B we show, that the above 
used definition of the free energy is accurate enough for all the p(k)s contributing to 
the Z significantly, as the uncertainty coming from the prescriptions for the JjS is 
0(l/vT), i-e. it disappears in the L — > 00 limit. 

2. The approximations in evaluating the Z . This group of questions is completely analo- 
gous to the questions emerged in connection to the free Fermi gas, and the answers are 
similar too: in an appropriate cutoff procedure involving the limits L — > 00, AA; — > 
and A — > 00 (Appendix D) Stirling's formula can be applied in its (2.6) form, q(p) can 
be taken at po and ph$ and also the limits of the £ integrals can be taken to infinity. 
(The major point in this is that the variable r oc £/\/LAk, i.e. only a ~ 1/V LAk 
neighborhood of the po and ph,o plays any role.) 

3. The application of the formulas (3.18) (the convergence of (3.20)). The condition for 
this is that all the eigenvalues of the matrix K are of modulus less than one. In the 
Appendix C we show that this is true for any T > 0, and we show also, that the T — > 
limit of the sum Y)„ -K, n exists too. 

As we have seen, our method is strictly established in the L — > 00 limit only, but in a 
less strict manner we may say, however, that the 0(1) corrections are correctly given also 
for finite but large enough size too: although in that case the A — > 00 limit can not be 
completed (as the LAkp(A) ^> 1 would not hold for very large A), the contribution of the 
states above A would be suppressed anyhow due to the large energy. For large but finite 
system, however, one has to be careful with the T — > limit. The quantity AS is an 0(1) 
correction to the thermodynamic potential: 



(3.22) 
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As for T — > F min — > _E (with E being the ground state energy) lim T ^ AS is an en- 
tropy. For infinite L this is the residual entropy which can be finite, but for large but finite 
L liniT^o A<S should be zero as the ground state is non degenerated. The resolution of 
this contradiction is that for finite L our calculation breaks down at temperatures where 
LAkpo t h(0) ~ 1, and below this temperatures the 0(1) corrections gradually disappear. For 
this reason the L — > oo and L» 1 cases should be distinguished! 



The case of open ends 



The system, just as in the previously discussed case, is described by the Hamiltonian (3.1), 
but the quantization condition is different: now the ring is not closed, and the particles are 
reflected on the ends. We suppose, that if a particle with wavenumber —k arrives at the 
end at x — 0, it is reflected to have a wavenumber k while the phase of the wavefunction is 
shifted by <po{k): 

e -ikx — > e ikx+iMk) ( x ~0). (3.23) 



Similarly, a particle of wavenumber k arriving at the end at x = L is reflected to have a 
wavenumber —k while the phase of the wavefunction is shifted by </?l(/c): 



e 



ikx > e -ikx+2ikL+i(p L (k) (x ~ L) . (3.24) 



(If the system is closed by infinitely high potential walls, <fo(k) = = n, but with other 

choices of <fo(k) and v ? l(^) different types of ends can be generated. In Appendix E we 
discus a case, when a surface potential having bound states closes the chain.) The Bethe 
Ansatz equations of such a system read: 

2Lk, = 2nJ t - (fo(ki) - cp L (h) ~ E ^tan" 1 {^^j +2^ ( [ h±!^ j . (3.35) 

Here the Jj numbers are integers, kj 7^ ±fcj for j ^ i, and none of the kjS equals zero. This 
means, the real kj are different positive numbers. Now we suppose, the <^o(k) and v 9 l(^) 
phases do not generate surface bound states (i.e. states with complex ks) and we discus 
the distribution of the real ks (the generalization for the case of surface bound states is 
straightforward, as it is seen in Appendix E). 

To proceed we split up the positive k axis into Ak intervals, and introduce the densities of 
the particles and holes (p(k) resp. in the usual manner. Now the integral equation 

connecting these quantities is 

p(k) + p h (k) = a(k) +J2K(k, k') p{k')Ak' 

k' 

roo 

= a(k)+ K(k,k')p(k')dk' (3.26) 
Jo 
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with 



{) 2n\ L dk L dk Lc 2 + (2k)\ 



and 



K (k, k') = — 



2c 



+ 



2c 



2n \c 2 + {k - k'Y c 2 + {k + k>) 



Now the number of states characterized by the same p{k) function is 

k)+p h (k 
LAkp(k) 



n <*wo) 

k=0 k^O 



Here k = refers to the interval beginning at the origin, and 



(LAk(p(k)+P h (k))\ 
\ LAkp(k) J ' 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



just as previously. The contribution of the k — interval differs from those of the others as 
the k = wavenumber (what is right on the edge of the interval) is not allowed. (In a 2Ak 
interval containing the origin in the middle there are 2LAk(p(k = 0) + ~p h (k = 0)) k values, 
out of this LAk(~p(k) +p h (k)) — 1/2 is positive, thus this is the number of choosable ks in 
the k — interval.) Applying Stirling's formula we arrive at 



n [p(k)\ = 



Ph(k) 



\ \p(k)+P h (k), 



n <"(?(*)) 

fc=0 k 



(3.31) 



The free energy functional is 

F \p{k)\ = E (Le(k)p(k)^k -T\nu{p{k))) ~ § In ?(fc ^ (fc) 



k=0 



Po(k) „-/3e(fc) 



(3.32) 



The minimization of the free energy leads to the condition 

= e~ p < k) , (3.33) 



Ph,o( k ) 

with e(k) given by (3.6) (just as in the PBC case), and it yields a minimal value 

F min = J± £ ]n (l + e ^( fe )) a{k)Ak + ^ In (l + e "^(°)) 
P 2 

= -=j ^ In (l + e-^W) a(A;)dA; + - In (l + e"^^) . (3.34) 
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Finally the evaluation of the functional integral leads to 

Z = e-^ndn+SS (3.35) 



with 



AS = E ^ ( 3 - 36 ) 



where 



/•OO /'OO 

XT = / • • • / dh-'-dknX (3.37) 
jo io 



PoOl) + Phfiih) ' Poifo) + Phflifo) Po( k n) + Phfiikn) 
/■oo /-oo ]_ ]_ X 

= io "70 dh ' ' ' dkn 1 + #M ^ 1 + #M • • • 1 + e*(*.) * (Av " ' 

Using the actual form of a(k) and K(k, k') it is easy to see, that 

-^min -^min + AF + 0Q + 0L, (3.38) 

with 

AF = I In (1 + e -*>) + | f -^*L_ In (1 + ,-**>) * , (3.39) 

We note here the following. The macroscopic part of the thermodynamic potential -F min is 
independent of the boundary condition as it should be. AF is a consequence of the openness 
of the chain, but the nature of the surfaces is reflected in O resp. 4>l only. These three 0(1) 
corrections directly connected to the surfaces are given already by the usual thermodynamic 
treatment. The contribution of the saddle point fluctuations depends on the boundary 
conditions in its structure. To see the difference between the AS valid for the PBC and AS 
applying for the open ends we introduce the function 

OO 1 

K(k, k') = ^2 -K n (k, k!) (3.41) 
n=i 71 



with 
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POO POO 

K n (k,k')= •••/ dJfcx • • • dAv,_i x (3.42) 

J —oo J —oo 

■M k M h , ^^ (fei'^) , , * , fc , ••• , 1, . K(k n . u k') ' 



Using the form of K(k, k') and K(k, k') it is not hard to see, that with this notation 

/OO , . roo , . /*oo , . 

K n (k,k)dk, and K=l K n (k,k)dk+ K n (k,-k)dk, (3.43) 

-oo JO JO 

thus 

/OO _ POO POO 

K(k,k)dk, while AS = / K(k,k)dk+ / K(k,-k)dk. (3.44) 

-co JO JO 

We note that the considerations presented in Appendix C concerning the convergence of AS 
hold for AS too, i.e. also AS together with its T — * limit exists. 

IV. GENERALIZATION 



The repulsive 5 Bose gas is the simplest BA system in the sense, that the particles have 
no internal structure and do not form bound states allowing to describe the system by 
one set of real parameters (the wavenumbers) . In most of the BA systems, however, for the 
thermodynamic description one has to introduce many (in most of the cases infinitely many) 
sets of rapidities. These can be of different type. For example in the case of Heisenberg 
chain these variables are the centers of the strings of different length [9], for the repulsive 5 
Fermi gas one set gives the wavenumbers, and the others are strings connected to the spin 
state of the system, while in case of the Hubbard model there are three type of rapidities: 
the real wavenumbers, the centers of the strings connected to the spins and the centers of the 
strings connected to the bound pairs [10]. (Here we take granted, that the string hypothesis 
works, as all in the known cases, where independent check is possible, it gives the correct 
result [11].) These systems are much more complicated than the 5 Bose gas, nevertheless it 
seems, that our calculation is generalizable for a larger class of them. Now we outline this 
procedure. 

For the sake of simplicity we denote all of the rapidities by k, this will cause no confusion. 
To each set of rapidities particle and hole densities (p^ n \k) resp. p^\k)) can be defined. 
Consider a system in which these satisfy a set of integral equations of the type: 



p in) 



{k)+pt\k) = a n (k) + J2 I K n:m (k,k')pW(k')dk'. (4.1) 



The number of states described by the same set of densities is now supposed to be 



16 



(4.2) 



and the free energy functional to be minimized is of the form 



p {n \k)] = L^J (e n {k)p {n \k) - Ts ( P {n \k))) dk 



(4.3) 



where e n (k) is the energy of an object of type n and rapidity k. The minimization leads to 
the condition 



p { o\k) 



= e 



-Pe n (k) 



with 



e n (k) = e n (k) - TJ2 /in (l + e~^^) K m>n (k', k)dk' . 



The minimal value of the free energy functional is 

F min = -TL W In (l + e~^) a n (k)dk . 



(4.4) 



(4.5) 



(4.6) 



The contribution of the states near to the one minimizing the free energy functional can be 
calculated by the functional integral method described in the previous section. Through a 
very straightforward calculation one finds, that the contribution of these states to the free 
energy given as the logarithm of the partition function is 



with 



where now 



AF = — TAS , 



as = e-^, 



(4.7) 



(4.8) 



^ n = E"--E / ••• / dh---dk n x 

mi m„ J 

1 1 

i + e fie mi (k 1 ) Km rrn 2 {k l , fc 2 ) _— g— . . . 



(4.9) 



1 



Kmn-rrn (k n , k\) . 



The above calculation is a formal generalization of the procedure applied for the 5 Bose gas, 
and it involves the same kinds of approximations too. For this in any case one should check, 
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if the conditions are met. Here we emphasize one point: Stirling's formula is applicable 
only if for all p^ a \k) giving significant contribution in the integral Lp^(k)Ak ^> 1 and 
Lp^\k)Ak ^> 1. This can be made true taking the L — > oo limit only in the case, if the 
Po\k) > and p]fy(k) > for all n and k. 

Our result has a meaning only if the sum defining AS converges. In any known system 
the kernel K mm /(k, k') is rather complicated, and checking the convergence may encounter 
difficulties. In some simple cases, however, this check is doable: For example in the case of 
the Heisenberg chain for T = the density of 1-strings (real rapidities) remains finite only, 
all other densities disappear, and the procedure of Appendix C can be applied. This gives 
the result that in the case of finite magnetic field the sum is convergent for the complete 
antiferromagnetic region (the anisotropy g > — 1 in the Hamiltonian (4.10)), but it is not 
convergent in the critical region (1 > g > — 1) for zero field. This renders it likely, that the 
sum is convergent for finite temperature too for finite field for any g, or even for zero field 
if g > 1, but does not support any guess for finite temperature and no field if 1 > g > —1. 

Even if the AS exists, to decide if (4.7) is correct for a given system may need further 
considerations. In deriving (4.7) we assumed that only (4.1) constrains the densities and 
that in principle all of the states are so described, i.e. (4.2) gives correctly the number of 
states described by the same density. These assumptions, however, may not be true even in 
the best known cases. Two kinds of problems may arise: 

1. In addition to (4.1) there are other constrains on the densities too. For example in 
systems of 1/2 spins the total number of the rapidities must not exceed the half of the 
number of sites/particles. (This expresses the fact, that the BA equations describe 
the states of total spin S z > only, the S z < states are obtained by reflection) 

2. Not all of the states are automatically described by the possible densities. For example 
for SU(2) systems the BA equations describe the states of highest weight only, the 
others are obtained from these by further manipulations. It is also possible, that in 
the system there are more vaqua, and more sets of excitations, and this multiplicity is 
not taken into account neither in the usual, nor in the present description. 

In the following we discuss these problems briefly in case of the Heisenberg chain. 

The XXZ Heisenberg chain of length N in magnetic field h is described by the Hamiltonian 

N 

Hxxz = E {SfSf +1 + SfSf +1 + gS*S* +1 - hS*) . (4.10) 

This model is Bethe Ansatz diagonalisable, but the BA equations give certain classes of 
the eigenstates only, and the others should be constructed by further manipulations. As 
the basic properties of the model are concerned four different case should be distinguished, 
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but from the thermodynamic point of view the following three antiferromagnetic cases are 
important. 

The — 1 < q < 1 planar case. In this case the BA equations give the S z = J2i Sf > states 
only: in all of the solutions the total number of turned down spins (spin-waves) is less than 
or equal to the half of the number of sites N/2. As it has been mentioned in point 1. above 
this imposes a constrain on the densities as they should satisfy the sum rule 

\>Y,hn)Jp {n \k)dk, (4.11) 

(n) 

with /(„) being the length of the string-type labeled by (n). Those densities not obeying this 
constrain are nonphysical, this can be manifested in negative /^(/^'s generated by (4.1). If 
in (4.11) the equality holds, than S z = 0, if the > sign is valid, than S z > 0, and the S z < 
states are constructed by reflection of the S z > ones. This means, that calculating the 
partition function the contributions of the S z > states have to be taken into account with 
weights 

1 + ex.p{-2(3S z h} , (4.12) 

i.e. all S z > states have additional contributions to the free energy 

-Tln(l + exp{-2[3S z h}) , (4.13) 

but no such contributions exist for the S z = states. For finite h the magnetization is 
macroscopic, S z oc L, thus if L — > oo, the above contribution to the free energy disappears, 
but for h = it remains finite —Tin 2. On the other hand in finite magnetic field all 
the densities contributing significantly to the functional integral obey the sum rule (the 
equilibrium densities satisfy (4.11) with the sign >), thus the correction (4.7) is certainly 
correct, provided it is convergent. This is also true for zero magnetic field, although the 
situation in that case is somewhat different. In zero field the equilibrium densities describe 
states with zero magnetization (i.e. they satisfy (4.11) with the sign =), and the functional 
integral also involves nonphysical densities (which violate the sum rule). Taking into account 
the symmetry of the functional integral it is clear, however, that the contribution of the 
nonphysical densities is the same as the contribution of the densities describing S z > 
(physical) states, i.e. the functional integral takes into account the S > states with a 
weight of two - just as it should be. 

The g — 1 isotropic chain. In this case in addition to the above constrain an other problem 
arises. The model has an SU(2) symmetry, and the BA equations describe the highest 
weight (S 2 = S Z (S Z + 1)) states only, the other states of the same spin length are obtained 
by means of the S~ operator. For this, if h ^ 0, the contribution to the partition function 
of one highest weight state should be weighted by 
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l-exp{-^(2g' + l)} f4 u) 

l-exp{-/3/i} 1 ' ' 

to get the contribution of the complete multiplet, i.e. to the free energy functional a term 

should be added. In the thermodynamic limit exp{— [3h{2S z + 1)} disappears, and we find 
an additional 0(1) correction 

Tln(l -exp{-ph}) (4.16) 

to the free energy. Note, that this is divergent if h — > indicating, that the correspond- 
ing correction in the h — case is of a different order of magnitude. (This divergence 
is a consequence of the thermodynamic limit, for finite N this term would behave as 
-Tln(2S z (h) + 1).) 

In the case of zero magnetic field all members of a spin multiplet are of the same energy, 
thus in order to get the contribution of a complete multiplet characterized by a spin length 
L (this L being equal to the S z of the highest weight member) a correction 

-Tln(2£+1) (4.17) 



should be added to the free energy functional. Although for the equilibrium density this 
correction is zero, we show in Appendix F, that it can be very large when deviating from 
the equilibrium density, thus its treatment requires a much more subtle procedure, what is 
beyond the grasp of the present work. We also give an estimation according to which this 
term is O(lniV), so we have to conclude that for the isotropic chain in zero field the leading 
correction to the macroscopic free energy may come from this term, not from the saddle 
point fluctuations. 

The g > 1 "easy axis" anisotropic chain. In this case - similarly to the g < 1 case - the 
solutions of the BA equations give all of the S z > states, and the S z < states are 
obtained by reflecting the S z > ones. In this sense the g > 1 and g < 1 are in close 
analogy and the conclusions concerning the 0(1) corrections made for g < 1 hold also for 
g > 1. Nevertheless in this case there is an additional problem. In this region the vacuum 
is twofold degenerated, this is manifested in the fact, that from the ground-state density 
of rapidities two different solution of the BA equations can be reconstructed. Also the low 
energy excited states can be grouped into two sets as being the excitations above one or the 
other vacuum [12], but it has not been studied yet, if such a degeneracy exists also in the 
thermodynamically important highly excited states. It is also not claerd yet, if our method 
takes this kind of degeneracy automatically into account, although this can be important, 
as it is expected to give an 0(1) correction too, which behaves like —Tin 2 as T — > 0. 
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V. SUMMARY 



In the present work based on the method deviced by Yang and Yang [3] we developed a 
functional integral method to calculate 0(1) corrections to the free energy of macroscopic 
BA integrable systems. In the Yang and Yang method the free energy of a system is written 
up as a functional of the momentum density, and this functional is minimized in order to 
find the actual value of the free energy. In terms of the grand canonical partition function 
the equilibrium density of momenta (at which the free energy is minimal) defines the states 
entering into the partition function with highest weights. The basic point of our calculation 
is that in evaluating the grand canonical partition function after the minimization of the 
free energy functional, (what actually gives the macroscopic part,) the contribution of the 
states near to the equilibrium (saddle point fluctuations) can be calculated by a Gaussian 
integral. To define this integral properly one needs to calculate the entropy entering into 
the free energy functional up to next to leading order. 

In addition to the technical problems the calculation of non macroscopic corrections to the 
macroscopic free energy rises some conceptional questions too. The Yang and Yang method 
has been developed to pick up the leading contribution only, thus in calculating further terms 
one has to see, that this refinement is meaningful, the method is accurate enough to calculate 
the next to leading contributions too. This involves two kinds of problems. The first is if it 
is possible at all to define an accurate enough free energy density in terms of the momentum 
density. (Questions of this type are discussed in Appendices A and B.) The other kind 
of problem is connected to the accuracy by which the macroscopic part itself is calculated: 
in systems, in which the number of available momenta is infinite, a cutoff procedure must 
be introduced even to calculate the macroscopic free energy. (The problem of this cutoff 
procedure and its resolution in the thermodynamic limit is discussed in Appendix D.) 

In order to avoid difficulties of special models we write up and discuss the method using the 
free Fermi gas. Our calculation reproduces the exact result, but we see that our method for 
the above mentioned implicit cutoff procedure is established in a strict mathematical sense 
in the thermodynamic limit only. 

Next the method is generalized for the repulsive 5 Bose gas with PBS. The structure of the 
BA equations for the rapidity densities of this system is (3.3) 




(5.1) 



with 



a(k) 



1 

27' 



K (k, k') 



1 



2c 



(5.2) 



2tt c 2 + {k - k') 



2 ' 



and the energy associated with a particle is 
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e(k) = k 2 -fi. (5.3) 

The 0(1) correction to the free energy we find is of the form —TAS with AS given in the 
form of an infinite sum (3.20-3.21) or equivalently in the form (3.41)-(3.44) 



3C 



AS= / K(k,k)dk, (5.4) 



where 



OO 1 

K(k, k') = ^2 -K n (k, k') (5.5) 

n=l 71 

with 

, . roc roc 

K n (k,k')= •••/ &i--4-iX (5.6) 

J —oo J — oo 



,/3e(fc') 

and the dressed energy e(fc) entering this formula being the usual one given by the equation 
(3.6) 

/OO , . 

K(k', k) In (l + e'^'A dk' . (5.7) 
-oo ^ ' 



We have shown also, that AS and its T — > limit exist and are finite. In the proof it is used 
that the particle density is a finite value fixed by the chemical potential /i. The "entropy" 
liniT^o AS is the residual entropy not equal to the entropy of the (unique) ground state 
what would be zero. This is a consequence of the fact that the T — > limit is taken after 
the L — > oo limit, and these limits do not commute. 

We also calculated the contribution of the saddle point fluctuations to the free energy of 
open end systems. We have found that it is slightly different in structure (3.44). It is —TAS 
with 

_ roc . roc . 

AS= / K(k,k)dk+ K(k,-k)dk. (5.8) 

JO JO 



In field theory systems of massive relativistic particles are important in which the particle 
number is not regulated by a chemical potential. A system of this type is described by the 
Lee- Yang model [13] in which e(k) = Mch(k) and the kernel K is negative. This system 
behaves somewhat differently than the 8 gas. There are strong arguments supporting the 
claim, that for PBC no 0(1) corrections shoud be present. Contrary to this our method 
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gives a AS of the type (5.4) for this case. For large T exp{— j3e(k)} is a constant for e(k) < T 
and is zero above. In this case AS oc — \n(T/M). If T — > 0, then e(&) — > e(/c), exp{/3e(/c)} 
diverges and A5 disappears. For the open end case we get a correction of the type (5.8) 
what also behaves differently as expected: due to the first term it diverges in the T — > oo 
limit. We note, that in a recent work [14] based on a different type of calculation it has been 
proposed, that the 0(1) correction in the open end case sholud be of the type of the second 
term in (5.8). 

One has to note here the following. It is possible to define some field theoretical models 
as certain limits of lattice models, for example the scaling limits of the Hubbard or the 
Heisenberg models are closely related to the SU(2) chiral invariant Gross-Neveu model 
[15-17]. Calculating the 0(1) corrections in the lattice models and taking the scaling limit 
afterward we expect lim^^o AS involve the residual entropy of the vacuum too. 

We also give a generalization of our result to other Bethe Ansatz systems. In a large class 
of models the densities satisfy equations of the type (4.1) 

P {n \k) +pt\k) = + E I K ntm {k,k')p {m \k')dk> . (5.9) 

m J 

For these systems we find, that the contribution of the saddle point fluctuations to the free 
energy is — TAS (4.7) with 

AS = Y,~K n , (5.10) 

n n 

where now 

•••/ dh---dk n x (5.11) 

I + e Pe mi {k 1 ) K ^n 2 {kl, fc 2 ) ^ + (fc 2 ) " " " J + {kn) ^m„ mi (K, h) . 

In this formula the energies e n (k) are connected to the e n (k) bare ones by the equations 

e n (k) = e n (k) - TJ2 I m (l + e~^^) K m , n {k\ k)dk' . (5.12) 

As this formula is a general one, its convergence should be checked in any special case. We 
also point out, that in the special cases additional problems requiring further considerations 
may arise, as we illustrate on the example of the XXZ Heisenberg chain. 
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APPENDIX A 



In this Appendix we want to examine the approximation (2.20) reading as 

E ex p I -0E (>I - } — > ex p I-' 9 ( ¥2 - + (ai) 

in more detail. Suppose that k is the mean value of the ks in the Ak interval, and that the 
Ak is small enough, so we may linearize around k. Thus to get the left hand side of the 
above formula we have to calculate 

{m \ 
2(7 1> > ( A2 ) 
i=i J 

where g = —(3kAk/N, the numbers are integers or half-integers (rij = (AT — l)/2 (mod 1)) 
satisfying -(AT - l)/2 < n x < . . . < m < n i+1 < . . . < n m < (N - l)/2 with N = 
L(p(k) + p (k))Ak and m = Lp(k)Ak ) and the J2{m} extends over all possible n; sets. This 
gives 

exp { -,3 (f - „) Lp(k)Ak} ft ^j^f 1>9> , (A3) 

what for small enough Ak yields 

-p (F - p) Lp{k)Ak + In u{p{k)) + ^ E ((^ - i + I) 2 - i 2 ) | • (A4) 

Evaluating the sum and inserting the value of N and m we arrive at 

exp [-p(k 2 -p- L P {k)Ak + \nu{p{k))) , (A5) 

It is not hard to see, that taking into account the quadratic nature of the spectrum would 
lead also to corrections not larger than 0((Ak) 2 ), i.e. the correction to the free energy density 
neglected in (2.20) ((Al)) is indeed small enough to disappear in the XX- • -)Ak — > / (. . .)dk 
limit. 



APPENDIX B 



Here we show, that the parity prescription for the parameters Jj (J; = (N + l)/2 (modi)) 
does not destroy the accuracy of the free energy functional. Let us consider two systems, 
one described by the system of equations 
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Lki = 2nJ, t - ^ 2 tan" 1 [^^j (Bl) 

with Jj being integer, the other with Jj being half-odd-integer independently from the par- 
ticle number. Let us distinguish between the parameters of the two systems by primes {k[, 
J-) and double primes (k", J"). (The real system is between this two: a primed solution 
is to be taken if the number of particles is odd, and a doubly primed one applies for an 
even number of particles.) Obviously the two kinds of solutions are in one-to-one correspon- 
dence: we consider a primed and a doubly primed solution one pair, if J" — J[ + \ for all 
%. Due to the 'Galilei invariance' of (Bl) the wavenumbers of the pairs are closely related: 
k'l = k\ + ir/L. For this we can describe the pairs by the same p(k) and Ph(k) if the Ah' and 
Ak" intervals are the same just shifted by n/L relative to each other. Let us denote the free 
energy associated to a p{k) in the two systems by F'[p(k)] and F"[p(k)], respectively. Their 
difference is 

AF[p(k)} = F"[p{k)\ - F'[p{k)\ ~ ]T 2kixp{k)Ak , (B2) 

k 

where on the right hand site we dropped the primes. For a general p(k) this can be of 0(1), 
but for those densities which contribute to Z it is much smaller: near to the equilibrium 
p(k) = po(k) +r(k), and as po is an even function of k, 

AF[p(k)} ~ ]T 2kir (r(k) - r{-k)) Ak . (B3) 

fc>0 

As in the functional integral the major contribution comes from the region r oc 1/y/L, for 
densities important in calculating the saddle point contributions AF[p(k)] oc l/VZ is a 
good estimation. This shows, that the difference in the free energies of the primed and the 
doubly primed systems disappears as L — > oo. As, however, (B3) is also an estimate for the 
error made if the parity of the numbers Jj is not chosen properly we may conclude that the 
prescription for the Js (Jj = (N — l)/2 (modi)) does not influence the 0(1) corrections. 



APPENDIX C 



In this Appendix we show, that all the eigenvalues k of the matrix K of (3.18) in the case 
of the 5 Bose gas have a modulus less than unity. To do this we use the formula 

maxlnkl = lim -ln|TrK n |. (CI) 
As TrK n is definitely positive for all n the eigenvalue of largest modulus is positive, i.e. 
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maxln \k\ = In /t max (C2) 

As we apply our formulas after the A/c — > / dfc limit is taken, we may use TrK™ = /C n , 
that is 

In /t max = lim — In K7 1 . (C3) 

For T > 0, using the relations 

/oo /-oo j 2uC 

-oo • " loo ^ ' ' ' d *<»-l>*(*' fel) • • • *(*<n-l), *0 = ^ {nc)2+{k _ kT , (C4) 

and 

/C ra is overestimated by the formula 

K n < / PQ(fc) V" 1 1 jV/g (C6) 



This yields 



«max < ( 7TT ~ 7TT J , (C7) 



what for T > is indeed less than one. This proves the applicability of the formulas (3.18) 
and the convergence of the series J2 n ( an d also proves, that the Neumann series of 
Eg. (3.7) converges). 

By a slight modification of the above estimations one can also prove, that the 

lim Y -KJ 1 (C8) 



also exists. For T = 

Po(k) f 1 if \k\ < fcj 



Po(k) + p h ,o( k ) ^° otherwise, 



(C9) 



where fop is a finite wavenumber, under which all, above which none of the states are filled. 
Observing, that 
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we see, that 



/Kp rltp 
■■• dk x --- dh(n-i)K(k, h) ■ ■ ■ K(k {n _ 1)} k) < 
-kp J —hp 

•••/ dk 1 ---dk (n _ 1) K(k,k 1 )---K(k (n . 1) ,k) = (CIO) 

-co J — oo TTTIC 



i-E^<-E^— • (en) 

t^o ^ n ire „ n 2 6c 



We note, that if c — > 0, i.e. if the kernel is a 5-function, both for T > and T — > our 
estimation blows up. 



APPENDIX D 

In this appendix we discuss the questions connected to the cutoff procedure involved in the 
evaluation of the free energy. First we notice, that this problem is rather a problem of 
the accuracy of the macroscopic part of the free energy. To see this consider the partition 
function Z of a system. The macroscopic part of the free energy is defined as 

i^mac = Lf , where f = -T lim - In Z . (Dl) 

The next to leading correction to the macroscopic part is 0(L°) if 

(AS =) hm In (ZeP F ) (D2) 

is finite but zero. If so, the free energy defined through the logarithm of the partition 
function behaves for large enough L as 

F = -T\nZ = Lf -TAS . (D3) 

In our case in order to apply Stirling's formula we have to introduce the cutoff A in the k 
space, and for the same reason we have to calculate sums instead of integrals. This way the 
partition function we obtain is of the form 

Z L (A, Ak) = e -/ 3i? min(L,A,Afc)+AS(A,AA : ) = e -/3L/ min (A,Afc)+A5(A,Afc) _ 



To make our reasoning simpler, for the time being we suppose, that it is accurate enough to 
replace the summations on AA; by integrals. This way we have 
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Z L (A) = e -/3i/ min (A)+AS(A) _ ( D5 ) 

Now 

-T lim y\nZ L = f min (A) (D6) 

from which the free energy density / is obtained through a next limit 

/ = lim / min (A) . (D7) 

A^oo 

This leads to difficulties in filtering out the next to leading corrections, as 

lim In (z L (A)e? Lf ) - (+or-)oo like L(f - / min (A)) , (D8) 

L— >oo V / 

i.e. / — /min(A) hides the correction we want to get. (We have to emphasize, that the value 
of the saddle point corrections is not effected by this, nevertheless we have to see the leading 
order term more accurately than the correction we expect.) 

To resolve this problem we propose the following. Taking larger and larger L allows taking 
larger and larger A, thus the two limits can be synchronized: a A(L) can bee chosen so, that 

A(L)— oo, (D9) 

while the condition for applying Stirling's formula within the cutoffs is met, i.e. 

LAA;p (A)>l. (D10) 

It seems plausible, that if it is possible 

• to choose A(L) so, that L(f — / min (A(L))) — > 0, 

• and take also Ak — ► so that the difference between the sums and integrals disappears, 

while (D10) holds, than the next to leading order correction to the free energy is 
liniA^oo AS (A) indeed. 

In the following we argue, that for the repulsive 5 Bose gas one can construct an appropriate 
cutoff procedure. First let us consider (D10). Due to the BA equations 

p -/fe(A) 

Po(A) = * 1 + e _ MA) (1 + 0(K(A, 0))) ~ ae-^ „ ae -Pe W ^ (Dn) 
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thus we require 

LAkae~ peW -> oo . (D12) 
(Here we used also, that for large A e(A) = e(A) + 0(K(A, 0)).) 

The effect of introducing integrals instead of the sums can be estimated by an Euler- 
Maclaurin type formula. We find the most significant part is 

~ L(Ak) 2 -^ae-^ A \ (D13) 

i.e. we need 

L(Ak) 2 ae~^A -> . (D14) 

(We do not give details here, just note, that making the error due to (2.20) (Appendix A) 
to disappear fast enough leads to the condition (D14) too.) 

Finally we had to estimate L(f — / m i n (A(L))) but for this we have to specify the procedure. 
One possibility is simply to omit all the modes outside the ±A interval (corresponding to 
taking their energy equal to oo), but in this scheme / min (A(L)) does not converge in A to / 
fast enough. A procedure providing a much faster convergence can be constructed realizing 
that the particles of high energy behave as free ones. In this scheme the free energy of the 
system is built up of two parts: the contribution of the modes within the ±A interval is 
calculated using Stirling's formula (just as in the bulk of the paper), while the contribution 
of the modes outside the cutoffs is approximated by the contribution of free particles of 
energy e(k) with density of states p(k) + Ph{k) given by the BA equations. In this scheme 
the minimization of the free energy leads to a dressed energy given by the equation 



e A (jfc) = e(k) — T j A K(k', k) In (l + e"^^) dk' (D15) 

f) 



~ T {l A + Ia) K{ ^' k) M ( X + e ^ e{k ' ] ) dk> > 



and the free energy density is 



/min(A) = - J A \n{l + e-^ A(k) )adk- + jf °°J In (l + e ~^) adk . (D16) 

The leading part of L(f — / min (A(L))) is of the order of 

~ L J™ (in (l + e -^ (fc) ) - In (l + e"^ k) )) adk , (D17) 
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what through some straightforward manipulations and approximations leads to the condition 

Le"^ A >^->0. (D18) 

It is not hard to see, that it is possible to define an L — > oo, A — > oo and AA; — * limit so, 
that all the conditions (D12)(D14) and (D18) are satisfied. 

We have to note, that the above reasoning concerning the existence of appropriate cutoff 
procedure works for the 8 Bose gas only, for other systems it has to be reformulated, but for 
certain models it is also possible, that there is no need for such a cutoff procedure. 



APPENDIX E 



In the present Appendix we discus a case of the 5 Bose gas with open ends, in which surface 
bound states can be present. The system is described by the Bethe Ansatz equations 

2Lk 3 = 2nJ 3 - Mkj) ~ Mkj) - £ \ Stan' 1 (h^tl) + 2 tan' 1 (h±Jl\ \ . ( E l) 




Now we suppose, ip L = n corresponding to an infinitely high wall closing the chain at L, but 
we take 

k 

ip (k) = 7T - tan -1 - . (E2) 

7 

Also this corresponds to an infinitely high wall, but this wall is preceded by an infinitely 
deep, but also infinitely narrow potential well. An appropriate tuning of the width and 
depth of the well leads to the above reflection phase shift. (The effect of such a potential in 
case of a 5 Fermi gas is discussed in [7], and a similar case of a Hubbard chain in [8]) This 
potential can always generate at least one surface bound state: it is not hard to see, that 
for any distribution of the real ks (we denote them by Latin ks) (El) has also imaginary 
solution too corresponding to a surface bound state at the end at x — 0. Denoting this by 
in we find 



K 



In the presence of such an imaginary wavenumber the real ones satisfy the equations 
2Lkj = 2nJj + tan" 1 h-t(k jlK )- £ {2 tan" 1 (^^j + 2 tan" 1 (^^j } , (E4) 
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with 



t(k, k) = {2 tan- 1 ( *— ^ + 2 tan" 1 ( *±*^ ) . (E5) 



If the potential is strong enough, i.e. 7 > c (El) has also solutions with more than one 
imaginary wavenumbers. These describe more than one particles bound to the x = end. 
For two imaginary wavenumbers ki i2 we find 

I 7 / ¥ (7 + + *f (7) 2 + *, 2 

(E6) 

In general, a solution with z/ imaginary wavenumbers iK a , a — 1, . . . v exists, if 7 — (u — l)c > 
0. The m a s have the form 

«i = 7 - 5 ± , /c a = re a _i -S a (a = 2,...u) (E7) 

with all 5 a being exponentially small in L. In these solutions the real k set satisfies the 
equations 

2Lkj = 2nJj + tan -1 ^ — tj,(%) — £ | 2 tan ^ (^T^) + 2 tan_1 (^T^) } ' ( E8 ) 
with 

V 

*„(*) = E «a) ■ (E9) 

a=l 

The thermodynamic treatment follows the procedure described in the bulk of the paper 
with the difference, that (p (k) + (p L (k) should be replaced by — 2 tan _1 (&;/7) + t u (k), and the 
calculation should be repeated for all possible v. Now the minimal free energy at a given v 
(taking also into account the direct energy contributions of the imaginary wavenumbers) is 

F v = F min + AF + AF U (E10) 
with F min being the bulk value, AF given by (3.39), and 

AFu=0= 2^i 7TF ln ( 1 + e k) ) dk > AF I/> o = AF +X;e(i« a ), (Ell) 
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with e{m) being the formal extension of (3.6) to complex k. As the saddle point contribution 
is independent of the state of the surface finally we arrive at 



/3F V +AS _ e -/3(F min +AF)+A5 



u=0 



i/=0 



-PAF V 



(E12) 



with n being the maximal possible value of v (that is always larger than or equal to one). 



APPENDIX F 



In this appendix we try to estimate the effect of the 

-Tm(2£ + 1) 



(Fl) 



term in the free energy of an isotropic Heisenberg chain in no magnetic field. As C is the 
magnetization of the highest weight member of the multiplet (which is the one described by 
the BA equations) 



N , \ 

C=--Nj2kn)T,P in) (k)Ak. 

A (n) 



(F2) 



(n) 

For the densities p given by the minimization of the leading part of the free energy 
functional this is zero, thus 



C = -Nj2l (n) J2r {n \k)Ak, 

(n) 



(F3) 



with 



r (")(jfe) =pW(k) -p { o\k). 



(F4) 



In analogy with (2.18) new variables are introduced 



\ 2 P r{k) P %{k) 



(F5) 



leading to 



C = ~ E l (n) E 

(n) 



2 AAA; ff (^KoW ^) {k) 



(F6) 
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As the application of Stirling's formula is correct if 

NAkpfr\k) > 1 and NAkp$(k) > 1 , (F7) 

-just as in the case of free Fermi or 5 Bose gas - one has to make sure through a cutoff proce- 
dure, that these inequalities hold. As the main contribution of the saddle point fluctuations 
come from the ~ 0(1) region we may conclude that for the thermodynamically 

important states £ can be large. 

The contribution of the (Fl) can be estimated as follows. Denoting the leading part of the 
free energy of the system as a function of the magnetization by F(S), the corrected free 
energy 

F(£) -Tm(2£ + 1) (F8) 



has to be minimized. This leads to 



(F9) 



where the susceptibility 



X 



d 2 F(S) 



OS 2 



oc N , 



(F10) 



5=0/ 



This way the minimal value is 



F(0)-Tln(2JT X + l) 



(Fll) 



indicating that the next to leading term is O(lniV), what is much larger than the 0(1) 
contribution of the saddle point fluctuations. 
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